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We analyze experimental data of nuclear structure function ratios F^^ / F^^ for obtaining optimum parton 
distribution functions (PDFs) in nuclei. Then, uncertainties of the nuclear PDFs are estimated by the Hessian 
method. Parametrization of nuclear parton distribution is investigated in the leading order of Qs. The parton 
distribution are provided at = IGeV^ with a number of parameters, which are determined by a analysis of 
the data on nuclear structure function. From the analysis, we propose parton distributions at — IGeV^ for 
nuclei from deuteron to heavy ones with a mass number A ~ 208. 



1. Introduction 

The structure function F2 for a bound nucleon, 
measured in deep inelastic scattering (DIS) of 
leptons, differs from that for a free nucleon. Un- 
polarized parton distribution in the nucleon are 
now well determined in the region from very small 
X to large x by using various experimental data. 
There are abundant data on electron and muon 
deep inelastic scattering. In addition, there are 
available data from neutrino reactions, Drell-Yan 
processes, W production, direct photon produc- 
tion, and others. 

In this paper we determine the unbounded par- 
ton distributions and structure function for free 
nucleon. For this peropus, we use the Bernstein 
average of moments, as we used in ill2,.3 4 . The 
polynomials are chosen so that the range of x for 
which the experimental values of F2 are not de- 
termined, make only a small contribution to the 
averages. These experimental Bernstein averages 
are then fitted in sec. 2, using CERN subroutine 

*e-mail: Atashbar@ipm.ir 

te-mail: Mirjalili@ipm.ir 

■t e-mail: Alinaghi.Khorramian@cern.ch 



MINUTE 0, to the QCD predictions for the 
corresponding linear combinations of moments. 
By having free parton distributions we can obtain 
the bounded parton distributions to determine 
nuclear structure function. In section 3 we calcu- 
late the weight function which takes into account 
the nuclear modification. Finally our nuclear 
parametrization studies are summarized in last 
section. 

2. Unbounded parton distribution in n and 

X space 

The evolutions of parton distribution func- 
tions in moment space are given by 

MuAn,Q^) - e-"^"^^ X MuJn,Ql) , 

MdAn,Q^) = e-''"^^ x Afd„(n,gg) , 

+ ((1 - a„)S(n, Ql) - f3.^g{n, Ql))e-''+' , 
g(n,Q2) ^ ((l-a„)g(n,Q2) 
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Here a„ and /3„ and associate anomalous di- 
mension are as following 
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the parameter s and b in above equations are de- 
fined as 
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If we back to Eq. we need to know all of the 
moments of parton at — Qq- 

In phenomenological investigations of structure 
functions, for a given value of Q^, only a limited 
number of experimental points, covering a par- 
tial range of values of x, are available. Therefore, 
one cannot directly determine the moments. A 



method devised to deal with this situation is to 
take averages of the structure function weighted 
by suitable polynomials. We can compare theo- 
retical predictions with experimental results for 
the Bernstein averages. After extracting some un- 
known parameters which exist in the functional 
form of input parton distribution in moment 
space, we are able to determine parton distribu- 
tions in a;-space. The parameterizations of the 
parton densities at input scale of = 1 GeV^ 
are determined as 
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We can control some parameter at input scale in 
above equations by these constrains 



Uy{x, Q^)dx = 2 



dv{x, Q^)dx = 1 



{xYj{x,QI) + xg{x,Ql))dx = I . (5) 



By taking the moments of the above equation and 
inserting them in the Eq. (^), we access to all of 
the parton distributions in moment space. Then, 
by using the inverse Mellin technique as 



/(a;,Q2) = - / dzlmle'^^x' 

Jo 



(6) 



we can obtain the parton distribution in a;-space 
for free proton. In above analysis, we chose the 
A = 0.22 GeV and c=l.l for non singlet part and 
c=2.1 for singlet part [7]- 
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Figure 1. Structure function for free proton in 
g2 = 9Gel/2 



3. Nuclear parton distributions and struc- 
ture functions 

Available nuclear data on the structure func- 
tion are taken in fixed-target experiments at 
this stage, and they are shown in Fig. 1 as a 
function of and x = Q^/(2 m i'), where i' is 
the transferred energy to the target, m is the nu- 
cleon mass, and is given by = —q^ with 
the virtual photon momentum q. The initial nu- 
clear parton distributions are provided at a fixed 
Q^{= Qg), and they are taken as [H] 

ffix,Ql)^w,ix,A,z)Mx,Ql) , (7) 

where ffix^Q^) is the parton distribution with 
type i in the nucleus A and fi{x, Qq) is the corre- 
sponding parton distribution in the nucleon. We 
call Wi{x, A, z) a weight function, which takes into 
account the nuclear modification. This functional 
form should be considered by: 

(1 ^) 

w,{x,A,Z) = 1 + + 

+ biX + Cix'^ + diX^ + Cix"^) . (8) 

Because the valence-quark distributions in a nu- 
cleus are much different from the ones in the 
proton, wc should be careful in defining the 



weight function Wi{x, A, Z). 

For nuclear A at Qg we have 

u^ix, Ql) = A, z fn^MD+^Nd^Ml) 

dtix, Ql) = A, z fd^MD+^NM^^Ql) ^ 

j:^{x,Ql)=w^ix,A,Z)J:{x,Ql) , 
g''{x,Ql)^Wg{x,A,Z)g{x,Ql) , (9) 

4. Summery 

In the theoretical calculations, the nuclei are 
assumed as: 

4iJe,^Li,9Be,12C,14Ar,27^^^40^„^56^g^63^j,^107^^^ 

"85n,i3iXe,i9^Au, and ^°^Pb 
The initial nuclear distributions are provided at 
Q"^ = I GeV^ with the parameters in Eq. (jH}. 
To obtain some unknown parameters which 
appeared in Eq. ©, we can use the ratios 
= F^{x, Q^)/F^{x, Q2) to calculate 

/ ryA,data jjA.theory\2 

y^=Y ~ - (10) 

j ^ j ' 

where the experimental error is given by consid- 
ering the systematic and statistical errors. Infor- 
mation about the used experimental data is given 
in Table I, where nuclear species, references, and 
data numbers are listed. 



Nucleus 


Experiment 


Reference 


No. of data 


He 


SLAC-E139 


lUJ 


21 




NMC-95 


Ha 


18 


C 


EMC-88 




9 




EMC-90 


cni 


13 




SLAC-E139 


in 


17 




NMC-95 


ini 


18 




FNAL-E665-95 


ini 


11 


Ca 


EMC-90 


nni 


13 




NMC-95 


ini 


18 




SLAC-E139 


in 


17 




FNAL_E665-95 


m 


11 



Table. I Nuclear species'references, and data 
numbers are listed for the used experimental 
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data with > 1 GeV'^ . 

With these preparations together with the 
CERN subroutine MINUIT V, the optimum pa- 
rameter set is obtained by minimizing ■ The ex- 
perimental data are taken from the pubhcations 
by the European Muon Collaboration (EMC)[n] 
at the European Organization for Nuclear 
Research (CERN), E139 Collaborations [IT] at 
the Stanford Linear Accelerator Center (SLAC), 
the New Muon Collaboration (NMC) ^ at 
CERN, and the E665 Collaboration [El at the 
Fermi National Accelerator Laboratory (FNAL). 
The used data are for the following nuclei: helium 
(He), carbon (C) and calcium (Ca). In Table IL 
we present the /d.o.f for helium, carbon and 
calcium as an example. 
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Nucleus 


No. of data 


X'id.o.f 


He 


39 


1.12 


C 


68 


1.49 


Ca 


59 


3.86 



Table. II The values, resulted from fitting, 
for different types of nucleus . 

In Fig. @ we present F^''(x,Q^)/F^{x,Q'^), 
F^ix, Q2)/Ff (x, Q2) and "(x, Q^)/F^{x, Q^) 
as a function of x and for = 5 GeV"^ . As it 
can be seen , these ratios are in good agreement 
to the available experimental data. 
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Figure 2. EMC effect for He, C, and Ca nucleus . 



3. S. Atashbar Tehrani, Ali N. Khorramian 
and A. Mirjalili, Commmi. Thcor. Phys. 
43(2005)1087. 

4. Ah N. Khorramian, A. Mirjalih, S. Atashbar 
Tehrani, Int. J. Mod. Phys. A 20(2005)1923. 

5. M. Cluck, R. M. Godbole, E. Reya, Z. Phys. 
C41(1989)667. 

6. F. James, CERN Program Library Long 
Writeup D506. 

7. M. Gliick, E. Reya, and A. Vogt, Z. Phys. C 
48 471 (1990) 

8. M. Hirai et ai, Phys. Rev. D64,(2001)034003; 
M. Hirai et a/., Phys. Rev. C70, (2004)044905 

9. EMC, J. Ashman et al, Phys. Lett.B 202, 
603 (1988) 

10. EMC, M. Arneodo et at, Nucl. Phys. B 333, 
1 (1990) 

11. SLAC-E139, J. Gomez et al, Phys. Rev. D 
49, 4348 (1994) 

12. NMC, P. Amaudruz et ai, Nucl.Phys.B 441,3 
(1995); NMC, M.Arneodo et al.,ibid. B 441, 
12 (1995) 

13. E665, M. R. Adams et ai, Z. Phys. C 67, 403 
(1995) 



